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Abstract 



The six-vertex model with domain wall boundary conditions is consid- 
ered. A Fredholm determinant representation for the partition function 
of the model is obtained. The kernel of the corrtesponding integral op- 
erator depends on Laguerre polynomials. 



1 Introduction 



In the present paper we consider the six-vertex model (ice model) on the two-dimensional square 
lattice with the special boundary conditions of "domain wall" ]IJ. First the six- vertex model was 
solved in the papers [J, where the partition function and bulk free energy were found for the 
case of periodically boundary conditions. An explicit expression for the partition function with 
domain wall boundary conditions on the lattice of finite volume was found by A. G. Izergin Q. 
Izergin's representation contains a determinant of N x N matrix (where N 2 is the number of 
vertices of the lattice). Thus, the problem of the thermodynamic limit of the partition function 
can be reduced to the study of the asymptotic behavior of the determinant at N — > oo. 

The goal of this paper is to obtain new representation for the partition function, appropriate 
for the asymptotic analysis. A Fredholm determinant representation seems to be suitable for this 
goal. Let us explain in brief the main idea of the method suggested. 

Let a matrix M. of the size N is given, whose entries Aijkif) depend on some parameter t. 
Determinants of this matrix det A4 generate a sequence of functions /at(£) for different N. As 
a rule the question on the properties of the limiting function /oo(t) = limTv^oo /jvCO is rather 
complicated. In spite of det.M is a finite sum of finite products, nevertheless such "explicit" 
expression for large N becomes extremely complicated for analysis. 

Suppose however that we can replace det Ai with a Fredholm determinant 

detM = det(/ + V) 
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of an integral operator, whose kernel parametrically depends on t and N: V(x,y) = V^{x,y\t). 
The Fredholm determinant is an infinite sum, and each term of the last one can be presented as 
a multiple integral, for example, 

det(7 + V) = exp ( ^2 J Vjv(xi, x 2 \t) ■ ■ ■ Vn( \t) d x . 

While at finite N the computations of the mentioned integrals are very complicated, at iV — ► oo 
the existence of large parameter in the integrands often allows to estimate not only each integral, 
but the series as a whole. Thus we have the possibility to obtain an estimate for the determinant. 
Hence, we can say that for iV — > oo representations in terms determinants of finite size matrices 
become "less explicit" , while Fredholm determinants turn to be "more explicit" . 

Transformations of finite determinants into Fredholm ones were found to be useful, for instance, 
in asymptotic analysis of correlation functions of quantum integrable models, asymptotic behavior 
of orthogonal polynomials, the spectrum of random matrices etc. We hope that in the case under 
consideration the Fredholm determinant representation will permit eventually to solve the problem 
of the thermodynamic limit of the partition function of the six-vertex model[[ 

The content of the paper is as follows. In the next section the necessary information on the 
six-vertex model and Izergin's formula for the partition function are given. In the third section we 
transform this formula to the Fredholm determinant representation. In the last section we discuss 
the outlook for application of the obtained formula to the asymptotic analysis. 



2 The six-vertex model 

Consider a statistical system on a square lattice of the size N x N, associating with each edge 
some variables (classical spin), taking values ±1. Usually these variables are denoted by arrows, 
identifying (f, — ►) with 1, and (J,, <— ) with —1. Thus, in general case there are 16 different 
configurations of arrows in each vertex. The six- vertex is specified by the condition that the number 
of arrows entering and leaving each vertex coincide. The types of possible vertices are shown on 
Fig. 1. With each configuration we associate a statistical weights. Generically the last ones can 
depend on the position of the vertex. We consider the model, which is symmetric with respect to 
simultaneous reversal of all arrows. Thus, we have 3iV 2 statistical weights Wjk = {ajk, bjk, Cjk}, 
where indices j, k = 0,1,..., N — 1 numerate vertices. The domain wall boundary conditions 
correspond to the entering arrows on the lower and upper boundaries, and leaving arrows on the 
left and right ones. 

1 Recently the thermodynamic limit of the six- vertex model was studied by use of Toda-chain equation || and 
matrix models methods (f| . The approach suggested in the present paper is completely different from the mentioned 
ones, and we hope that it is also worth attention. 
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bjk 
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Figure 1: Verteces of the six- vertex model and their statistical weights. 



The partition function is given by 



N-l 

Z N = II W 3k- 

{C} j,k=0 



(2-1) 



Here the sum is taken with respect to all possible configurations C on the whole lattice. 

The six-vertex model is exactly solvable, if 3iV 2 weights bjk, Cjk} are parameterized by 

2N + 1 variables Xj, (j = 0, 1, . . . , N - 1) and 77: 



a jk = <p(\j - €k + v), bjk = tp(Xj - Cj k = ip{j}). 



(2.2) 



Parameters Aj are associated with horizontal lines, £ 3 - — with vertical ones (see. Fig. 2), 77 is an 
arbitrary complex constant. Function (p(x) is defined up to normalization factor and is equal to 



(p(x) 



sinh(x) 



or 



ip(x) 



x. 



(2.3) 



In the first case the model is equivalent to the quantum XXZ magnet, in the second — to the 
XXX magnet. 

The Izergin's formula for the partition function on the lattice of finite volume with the domain 
wall boundary conditions has the form 



N-l 



II <p(K - €b + r])ip(X a - £ b ) 
Z N = = ^ — — det 



<p(v) 



I! ip(X a - X b )tp(£, b - £ a ) \<p(\j - Ck)^f(Xj - Ck + v) 

N>a>b>0 



(2.4) 



We would like to draw attention of the reader that here and further we numerate indices of matrices 
starting with zero, but not with one: j, k — 0, 1, . . . , N — 1. 
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Figure 2: Parametrization of the edges and verteces of the lattice. 



We are interesting in homogeneous limit, when the statistical weights do not depend on the 
position of the vertex. In other words Xj = A, £_j = £, j — 0, 1 . . . , N — 1. The corresponding 
limit in (12. 41) can be obtained via 



A(A)A({) 



i <9 j ' a* 



■$(A,0 



(2.5) 



Here $(A,£) is some enough smooth two- variable function, A (A) and A(£) are Van der Monde 
determinants of variables {A} and {£} respectively. 

In the rational case, when (p(x) = x the direct use of (|2.5|) gives 



7 {XXX) 
J N 



N 2 



N-l 

I! (m!) 2 

m=0 



■det 



(j + k)\(i 



-j-k-i 



> + v) 



-j-k-i 



(2.6) 



where v = X — £. 

In the trigonometric case it convenient first to use variables Uj = e 2Aj , Vj = e 2 ^ and q = e v . 



Then ( [2.4|) takes the form 



N-l 



7 {XXZ) 



N-l 

n 

a=0 



I! sinh(A a - £ 6 + r}){u a - v b ) 



n (w a - u b )(v b - v a ) 

N>a>b>0 



det 



1 



Uj - w fc quj - q L v kj 



(2.7) 



Now the equation ( |2.5| ) can be easily applied, and after simple algebra we obtain 



y(XXZ) _ Nv 

N — e 



sinh(zy) sinh(z/ + rj) 



N- 



N-l 

II (m!) 2 

m=0 



■det 



{j + k)\ qi- k+1 (j + k)\ 

sinh^+V) ~ sinh J+fc+1 (z/ + ?7) 



(2.8) 
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Here as in fl2.6| ) v = A — £, and q = e v . 

Below we shall focus on the trigonometric case (|2.8|) , as the most general one. As usual, the 
answer for the partition function Z^ xxx ^ can be obtained in the limit v = ei/, rj = er], e — ► 0. 

Our goal is to transform (|2.8| ) to the expression, suitable for the asymptotic analysis at N — > oo. 

There are factors with trivial asymptotics in the equation ( |2.8| ): e and sinh(z/) sinh (v+rj) . It 
is also easy to evaluate the behavior of the factor ]lm=o( m 2 at iV — > oo. The main problem is the 
determinant. In the next section we consider some special transforms of (|2.8|) , which result finally 
to the new representation for the partition function, containing the determinant of an integral 
operator. 



3 Fredholm determinant 

As the first step, it is convenient to extract, for example, siT\h~^~ k ~ 1 {u) out off the determinant 
(gg ). Then (fZ|) becomes 

z^ xz) = e N ^ N y+v) det[A _ qtQTATQ -^ ( 3. 1} 

11 (m!) 2 

m=0 

Here 

sinh ( v ) 

A jk = (j + k)\, T jk = 6 jk lP, Q jk = S jk q 1 , t= , Y (3.2) 

smh(^ + 77) 

Note. Similarly we could extract terms sinh^^Yz/ + v)- Then, up to the replacements 
t — > i -1 , g — ► and trivial common factor, the arising expression coincides with ( |3 . 1| ) . 
Now we can transform (13. II) as follows: 



det[A - qtQTATQ- 1 } = det A ■ det[J - U], where C7 = qtA~ l QTATQ~ l . (3.3) 

The advantage of this method is that det A cancels the product of factorials in the denominator, 
since 

N 

det A= n( m! ) 2 > ( 3 - 4 ) 

(see Appendix 1). Thus our problem reduces to the derivation of the inverse matrix A' 1 and 
calculation of the product U = qtA~ 1 QTATQ~ 1 . The inverse matrix A' 1 can be easily written in 
terms of Laguerre polynomials, therefore below we recall their basic properties. 
The Laguerre polynomials L n (x) are defined by 

L n (x) = -^{x n e~ x ) (3.5) 
n! dx n 

and generate the orthogonal system on the half- axis R + with the weight e~ x : 

L n (x)L m (x)e~ x dx = 5 nm . (3.6) 
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The key object, determining the matrix A 1 , is the kernel 



K n (x,y) = J2 L k(x)L k (y). 

k=0 

One more useful representation for K n (x,y) follows from Christoffel-Darboux formula: 

n + 1 



(3.7) 



K n (x,y) 



x-y 



(3.8) 



L n+1 (x)L n (y) - L n+1 (y)L n (x) 

Obviously, for arbitrary polynomial 7r m (x) of degree m, less or equal to n, holds 

■K m (x)K n (x,y)e~ x dx = ir m (y), 0<m<n. (3.9) 



To prove (|3.9j), it is enough to substitute to the l.h.s. the expansion of vr m (x) with respect to 
Laguerre polynomials 7r m (x) = J2k"=o c kLk(%), and to use the definition (p.7|). 

Being a polynomial of variables x and y, the kernel K n (x, y) can be presented in the form 



it 

K n (x,y)= J2 K%y y k , 



where 



K 



(n) 



j,fc=0 

1 &> +k 



K n {x,y) 



(3.10) 



(3.11) 



jk ' ]\k\dxWy k ' 

Now we are in position to formulate the theorem on the inverse matrix. 

Theorem 1 Let Ajk = (j + k)\, j, k = 0, 1 . . . , iV — 1. Then the inverse matrix has the entries 

(A- 1 )* = Kj;- 1 * . (3.12) 
Proof. Using the integral representation for (j + £;)!, we have 

I Qk N-l ^ Ql roo 



N-l 

2 / AjiK lk 



.(JV-l) 



1 u f°° 



(3.13) 



x=y=0 



The sum with respect to / is the Taylor series (in fact, since K^_i{x, y) is (TV— l)-degree polynomial 
of x, the mentioned series turns to the finite sum). Hence, 



N-l 



i a* 



roo 

/ s J e~ s K N _ 1 (s,y)ds 
Jo 



i=o kl °y k 

Due to (|3~9|) the integral in (|3.14j ) is equal to yK Thus, we arrive at 



(3.14) 



y=0 



N-l i ofe 

V A-iK^-V - - — 

to 3 ~ k] - 9y k 



yj 



Ijk- 



(3.15) 



y=0 



The theorem is proved. 

Note. For generalization of the relationship between Hankel matrices and Christoffel-Darboux 
kernels see J7j. 

The obtained explicit expression for A" 1 allows us to compute the matrix U Q3.3|) . Indeed, 



U 



1 & 
k\ dy k 



k N-l 



E^V-^ 1 — jf s^e-K N . l{x , y) ds 



(3.16) 



x=y=0 



Similarly to the proof of the Theorem the sum with respect to / easily can be computed and it 
gives 

qt d k roo 



Ujk = \\-QkJ xj ^ q XKN -^ xqt ^ y"> dx 



(3.17) 



Thus, the partition function Z^ X turns to be proportional to the determinant of the matrix 
5jk — Ujk ( |3.17| ). To compute the last one it is convenient to use the equation 



det(J -U) = exp{tr log(J - U)} = exp I - - tr U n \ . 



n=l 



n 



Hence, we need to find the traces of powers of the matrix U. It is easy to see that 
(qt) n d k [°° 



k\ dy k Jo 



POO sr-^n 

/ x> l (q- l t) k K N _ 1 (x 1 qt, x 2 q~H) ■ ■ ■ K N ^(x n qt, y)e~^m= 1 d T 
Jo 



X 



(3.18) 



(3.19) 



y=0 



Indeed, for n = 1 ( |3.19| ) coincides with ( 3.17 ). Assuming that ( |3.19| ) holds for U n , we immediately 
obtain 



JV-l 

E (u n hUi k 

1=0 



(qt) n+1 d k f°° 



k\ dy 



POO 

: / x J 1 (q-H) k K N - 1 (x 1 qt,x 2 q~ 1 t) 
• Jo 



■■■K N ^ l (x n qt } x n+1 q~H)K N _ 1 (x n+1 qt } y)e-^^i Xm d n+1 x . (3.20) 

y=0 

(Here the sum with respect to I again turns to the Taylor series). Hence, 

f°° i 1 v^™ 

ti(U n ) = {qt) n / K N _ l (x l qt,x 2 q-H)---K N _ l (x n qt,x l q- l t)e-^m=^d n x. (3.21) 
Jo 

The obtained expression ( |3.21| ) for the trace of U n allows to replace det (£,,•& — Ujk) with Fredholm 
determinant. In order to do this, we introduce integral operator I —V, acting on R + 



with the kernel 



[(I-V)f](x) = f(x)- V{x,y)f{y)dy, 

Jo 



N-l 



(3.22) 



V(x, y) = qtK N -i{xqt, yq-^e'^' 2 = qt £ L k (xqt)L k (yq-H)e~^/ 2 . (3.23) 

fc=0 



Due to 



V(x, y) also can be presented in the form 



V(x,y) = Q _ - t (L N ^(xqt)L N (yq-H) - L N {xqt)L N _ x (yq-H)) e~^ 2 . 
qx q y 



_qN 

Then due to ( ggg ) 
and we finally obtain 



trU n = trV n , 
det(5 jk -U jk ) = det(I-V). 



(3.24) 

(3.25) 
(3.26) 



4 Discussions 

Thus, we have found new representation for the partition function of the six-vertex model, con- 
taining the Fredholm determinant of the integral operator with the kernel ( |3.23|) , (|3.24 ): 
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J N 



N 2 



sinh(^ + rj) det(I-V). 



(4.1) 



Since the kernel V is degenerated, we always can turn back to the determinant of the finite size 
matrix 



N-l 



det II - qt J2 L k (xqt)L k (yq H). 



,-(x+y)/2 



det S 



k=0 



qi J Lj(xqt)L k (xq 1 t)e x dx) . (4.2) 



oo the Fredholm determinant representation seems 



However, for the asymptotic analysis at N 
to be preferable due to several reasons. 

First, the equation (|3.9| ) shows that the kernel K N (x,y) acts as identity operator on the sub- 
space of the polynomials of iV-th degree. Then one can think that for N going to infinity, Kn(x, y) 
goes to identity operator (delta- function). Of course, this limiting procedure requires a serious 
analysis of convergency. It is remarkable however, that for certain values of q and t such a naive 
interpretation of the kernel ( |3.23|) does holds, and the Fredholm determinant can be evaluated 
explicitly in the limit N — > oo. 

Second, similar integral operators with the kernels, depending on Laguerre polynomials, often 
appear in the theories of random matrices and random permutations || || . Some methods of the 
analysis of such operators can be directly applied to the case under consideration. Besides, at 
large iV the Laguerre polynomials can be approximated by Bessel functions, what also leads us to 
the integral operator, analogous to ones, considered in the papers |Tt], |TT], 

J N 

the form 



Third, in particular case of the partition function Z^ xx \ when q = 1, the kernel V(x, y) takes 



V(x,y) 



N 



x-y 



{L N - X (xb)L N (yt) - L^xt^^iyt^e-^' 2 . 



(4.3) 



The integral operator with the kernel ( |4.3|) belongs to the class of integrable integral operators 
p~3f . The corresponding Fredholm determinant turn to be r-function of classical exactly solvable 
equation and can be evaluated via the matrix Riemann-Hilbert problem methods. 



We are planning to present the detailed asymptotic analysis of the partition function of the 
six-vertex model with domain wall boundary conditions in forthcoming publications. 



The author thanks A. R. Its for numerous and useful discussions. The work was supported in 
parts by RFBR, Grant 99-01-00151 and INTAS-99-1782. 



A Calculation of the determinant 

One of the method to proof (|3.4j) is to use well known formula for Cauchy determinant 

I! (a a - a b )(p b - Pa) 



N>a>b>0 



1 



a,b=0 

Here {a} and {/?} are arbitrary complex. Dividing both sides of (|A.1|) by Van der Monde deter- 
minants A(a) and A(/3) and taking homogeneous limit atj — > a, flj — > /3, j = 0,1, . . . , N — 1, 
we obtain due to ( |2.5| ) 



m=0 



what immediately implies 

Another method to compute det A is based on the application of Laguerre polynomials. Below 
we present this method, since it allows to find A' 1 as well as to prove 

Using the integral representation for factorial, we have 



det A = det (j e~ s s j+k dsj . (A.3) 
Denote integration variable in the j-th row as Sj. Then we come to the multiple integral 

POO C CXj 

detA= ds - ■ -dsN^det (e~ Sj 'sj +fc ) = / ds ■ ■ ■ ds N -i ]J (e~ Sj s J ) A(s), (A.4) 

JO JO j=0 

where A(s) is Van der Monde determinant. Now we can use well known identity 

A(s) = det[p fc (s i )], (A.5) 

where Pk{s) is a system of arbitrary polynomials with the highest coefficient equal to 1: 

p k (s) = s k + o(s). (A.6) 

In particular we can choose the Laguerre polynomials 

p k (s) = L k (s) = (-l)V-^-( S fc e- s ). (A.7) 



The normalization in ( |A.7| ) differs from one in ( |3.5p by factor (— l) k k\ in order to provide (A. 6). 



Therefore we denote the polynomials ( |A.7| ) as L^. The orthogonality condition then takes the form 

roo _ _ 

/ L m (s)L„(s)e- s ds = (n!) 2 5 nm . (A.8) 
Jo 

Thus, the equation (|A.4|) becomes 

roo N - 1 

detA= ds ---ds N -i ]J (e~ Sj s j ) det[L k ( Sj )}. (A.9) 

Integrating inside the determinant, we obtain 

det A = det ( r e- s s j L k (s) ds) . (A. 10) 



Due to orthogonality of the polynomials Lk(s), the entries in ( [A. 10 ) vanish, if j < k. Thus we 
deal with the determinant of triangular matrix and, hence, 

N— 1 / poo „ \ 

detA=n(/ e- s s m L m (s)ds . (A.ll) 

m=0 ' 

It remains to observe that expansion of s m with respect to L m (s) has the form 

s m = L m (s) + polynomials of lower degree 

therefore due to ( |A.8|) 



N-l 

detA= I]( m! ) 2 , (A.12) 

m=0 

which ends the proof. 

The described method allows to compute not only det A, but the minors of the size — 1 as 
well. This provides the possibility to find A -1 . 
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